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We examined a model of partially screened gap region 
above the polar cap, in which the electron-positron plasma 
generated by sparking discharges coexists with thermionic 
flow ejected by the bombardment of the surface beneath 
these sparks. Our special interest was the polar cap heat¬ 
ing rate and the subpulse drifting rate, both phenomena 
presumably associated with sparks operating at the po¬ 
lar cap. We investigated correlation between the intrinsic 
drift rate and polar cap heating rate and found that they 
are coupled to each other in such a way that the thermal 
X-ray luminosity L x from heated polar cap depends only 
on the observational tertiary subpulse drift periodicity P 3 
(polar cap carousel time). Within our model of partially 
screened gap we derived the simple formula relating L x 
and P 3 , and showed that it holds for PSRs B0943+10 and 
B1133+16, which are the only two pulsars in which both 
L x and P3 are presently known. 
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Abstract. 


1. Introduction 


Revealing the problem of the drifting subpulses phe¬ 
nomenon seems to be a keystone for the understanding 
of the pulsar radio emission mec hanism. In the model 
of lRuderman .^TSutherlanj ( 19751 RS75 henceforth) the 
spark-associated subbeams of subpulse emission circulate 
around the magnetic axis due to E x B drift of spark 
plasma filaments. This model is widely regarded as the 
most natural and plausible explanation of drifting sub¬ 
pulse phenomenon, at least qualitatively. Despite its popu¬ 
larity the RS75 model is known to suffer from the so-called 
bindi ng energy problem (for review see lllsov fc Melrosd 
11996 1. In fact, the cohesive energy of surface iron ions 
were largely overestimated and the ’’vacu um gap” en- 
visioned by RS75 was impossible to form. iGil fc Mitral 
il 2 nmh and El fc Melikidza (1299211 revisited the binding 
energy problem and concluded that the formation of RS75 
gap was, in principle, possible, although it required a 
very strong surface magnetic field, much stronger than 
the dipolar component inferred from t he o b serv ed spin- 
down rate. Recent calculations of iMedin fc Lail ( 2006 1 
showed that iron chains are strongly bound in mag¬ 
netic field close to 10 14 G. This can only occur in ac¬ 
tual pulsars if their surface field is dominated by strong 
non-dipolar components. Growing evidence of such non- 
dipolar surface anomalies accumulates in the literature, 
both observationa l an d theoretical (for short review see 
e.g. llJrpin fc Gill . 12004 . Therefore, it seems natural to 
assume that the magnetic field at the stellar surface is 
dominated by the crust anchored small scale (sun-spot 
like) components, which can significantly alter the physi¬ 
cal condi tions for plasma product ion in the polar cap re¬ 
gion (see lGil. Melikidze. fc Mitr3 . 129921 for details). 


Even if the RS75 gap was possible to form, it would 
not automatically solve the mystery of drifting sub¬ 
pulses, since it predic ts to o fast a drifting rate (e.g. 
Deshoande fc Rankinl Il999ll. M otivated by this issue 


Gil. Melikidze fc Geppertl d299.lL GMG03 henceforth) de¬ 


veloped further the idea of the inner acceleration region 
(IAR henceforth) above the polar cap by including the 
partial screening due to thermionic ions flow from the sur¬ 
face heated by sparks. We will call this kind of IAR the 
’’partially screened gap” (PSG henceforth). The slow drift 
corresponds to E x B drift of spark generated electron- 
positron pairs, until the total charge density within the 
PSG reaches the co-rotational value. Since the PSG po¬ 
tential drop is much lower than in the RS75 model, the 
intrinsic drift rate is compatible with observations (for de¬ 
tails see GMG93). 

The PSG model allows a relatively high heating rate 
of the polar cap surface, compatible with observations, 
in contrast to RS75 gap, which produces too much heat. 
On the other han d , the alternative Space Charge Limited 
Flow model dArons fc Scharlemanill97 9tl predicts a much 
lower heating rate (e.g. Zhang fc Hardinet l200Cil . Thus, 
measuring the thermal X-ray luminosity from heated po¬ 
lar caps can potentially reveal the nature of IAR. This 
can help us to understand a mechanism of drifting sub¬ 
pulses, which appears to be a common p henomenon in 
radio pulsars. Recently IWeltevrede. Edwards fc Stapperd 
(299(1 WES06 henceforth) presented results of a system¬ 
atic, unbiased search and found that the fraction of pul¬ 
sars showing drifting subpulses is at least 55 %. They con¬ 
cluded that the conditions for drifting mechanism to work 
cannot be very different from the emission mechanism of 
radio pulsars. 

In or der to test diff erent avai la ble m odels of IAR in 
pulsars iZhang. Sanwal fc Pavlov! ( 120051 . ZSP95 hence- 
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forth) observed the best studied drifting subpulse radio 
pulsar PSR B0943+10 with the XMM-Newton. Their ob¬ 
servations were consistent with PSG formed in strong, 
non-dipolar magnetic field j ust above the surface of very 
small polar cap. Recently. Ilvargaltsev. Pavlov, & Garmirel 
(. 2006 . KPG06 henceforth) observed the X-ray emission 
from the nearby PSR B1133+16 and found that this case 
is also consistent with the thermal radiation from a small 
hot spot (again much smaller than the canonical polar 
cap). PSR B1133+16 is almost a twin of PSR B0943+10 
in terms of P and P values and, interestingly, both these 
pulsars have very similar X-ray signatures, in agreement 
with our PSG model (see Table 1). 

In this paper we generalized the treatment of ZSP05 
and developed a model for thermal X-ray emission from 
radio drifting pulsars. The model matches the observa¬ 
tions of PSRs B0943+10 and B1133+16, both in radio 
and X-rays. There is a number of pulsars with measured 
thermal X-ray radiation from small hot polar cap but un¬ 
fortunately their drifting properties are not yet known. 
This is likely to change in the near future due to new so¬ 
phisticated methods for analysis of intensity fluctuations 
in weak pulsars being developed (e.g. WRS06 and refer¬ 
ences therein). 


2. Charge depleted inner acceleration region 


The IAR above the polar cap results from the devia¬ 
tion of a local charge density o from th e co-rotational 
charge density 1 Goldreich fc JulianL Il969t GJ69 hence¬ 
forth) pcj = — ri ■ B s /27tc « B s /cP. As already men¬ 
tioned, growing evidence appears, both observationally 
and theoretically, that the actual surface magnetic field 
B s is highly non-dipolar. One can introduc e a scaling fac- 
tor b = B s /Bd to describe its magnitude dGil fe Sendvkl 
200 ft GS00 henceforth), where the enhancement coeffi¬ 
cient b > 1 and B d = 2 x 10 12 (PP_is) 1//2 G is the canon¬ 
ical, star centered dipolar magnetic field, P is the pulsar 
period and P_i 5 = P/ 10~ 15 is the period derivative. 

The polar cap is defined as the locus of magnetic 
field lines that penetrate the so-called light cylinder. 
Conventionally, the polar cap radius r p 


10 4 P 0 5 cm, and its surface area A nr ~ 2 x 10 8 P 1 "™ 2 


= 1.45 x 
^ cm 2 . 

In the case of non-dipolar surface field the polar cap area 
must shrink due to flu x conservation of t he op en field lines 
A pc B d = A p B s (e.g. ICheng fc Zhand. Il999l) . Thus, the 

~ 1 A„ n , regard- 


surface area of the actual polar cap A p = b j-l pc 
less of the actual shape of the polar cap. Consequently, one 
can write the polar cap radius in the form r p = b~ 0 5 r pc . 
However, one should realize that this is only a characteris¬ 
tic dimension of the actual polar cap. Indeed, the presence 
of strong surface magnetic field anomalies should make the 
shape of the polar cap quite irregular. 

For isolated neutron stars one might expect the surface 
to consist mainly of iron formed at the neutron star’s birth 
( Lai 12001 ; i Usov fc Melrosei . Il996 !. Therefore, the charge 
depletion above the polar cap can result from bounding 
of the positive I^Fe ions (at least partially) in the neu¬ 


tron star surface. If this is possible, then positive charges 
cannot be supplied at the rate that would compensate the 
inertial outflow through the light cylinder. As a result, 
a significant part of the unipolar potential drop (GJ69, 
RS75) develops above the polar cap, which can accelerate 
charged particles to relativistic energies and power the 
pulsar radiation mechanism. 

The cascading production of electron-positron plasma 
is crucial for l imitat ion of growing gap potential drop (e.g. 
R.S75. fGheng fc R,udermarill98(i CR80 henceforth). The 
accelerated positrons will leave the acceleration region, 
while the electrons will bombard the polar cap surface, 
causing a thermal ejection of ions. The ions will cause 
partial screening of the potential drop, which can be de¬ 
scribed as AV = r](2ir/cP)B s h 2 , where h is the height of 
the acceleration region, r] = 1 — Pi/P gj (CR80) is a shield¬ 
ing factor and p* is charge density of ejected ions. The 
gap potential drop is completely screened when the total 
charge density p = pt + p+ reaches the co-rotational value 
pGj = en.Gj, where uqj = 1.4 x 10 11 6P° 1 5 5 P _0 - 5 cm -3 is 
the co-rotational charge number density . 

GMG03 argued that the actual potential drop AV 
should be thermostatically regulated and the quasi¬ 
equilibrium state should be established, in which heating 
due to electron bombardment is balanced by cooling due 
to thermal radiation. The quasi-equilibrium condition is 
Qcool = Qheat: where Q coo i = crTf is a cooling power 
surface density by thermal radiation from the polar cap 
surface and Qheat = 7 m e c 3 n is heating power surface den¬ 
sity due to back-flow bombardment, 7 = eA V/m e c 2 is 
the Lorentz factor, n = hqj — n t = r/ncj is the number 
density of back-flowing plasma particles depositing their 
kinetic energy at the polar cap surface, rii is the charge 
number density of thermionic ions. It is straightforward 
to obtain an expression for the quasi-equilibrium surface 
temperature in the form 

T s = (6.2 x 10 4 K)(P_i 5 /P) 1 / 4 ? 7 1 / 2 6 1 / 2 /i 1/2 . (1) 


Let us now interrelate the accelerating potential drop 
AV and the perpendicular (with respect of the magnetic 
field lines) electric field A E which causes E x B drift. 
Following of RS75 we can argue that the tangent elec¬ 
tric field is strong only at the polar cap boundary where 
A E = 0.5AV/h = rj(jr/cP)B s h (see Appendix A in 
GMG03 for details). Due to the E x B drift the discharge 
plasma performs a slow circumferential motion with ve¬ 
locity v d = cAE/B s = rjirh/P. The time interval to 
make one full revolution around the polar cap boundary 
is P 3 ~ 2irr p /vd- One then has 


h _ 


P 2rjh 


If the plasma above the polar cap is fragmented into fil¬ 
aments (sparks) which determine the intensity structure 
of the instantaneous pulsar radio beam, then in princi¬ 
ple, the tertiary periodicity P 3 can be measured/estimated 
from the pattern of the observed drifting subpulses 
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Table 1 . Comparison of observed and predicted parameters of thermal emission from hot polar caps 


Name 

P3/P 

P3/P 

Lx/P 

x 10 3 

b 

rr i(obs) 

J- S 

Ti(pred) 

± s 

B d 

B s 

PSR B 

Obs. 

Obs. 

Pred. 

Obs. 

Pred. 

-^pc/-^bol 

O 

O 

t-H 

10 b K 

10 i2 G 

TcPg 

0943 + 10 

1.86 

37.4 

361“ 

n 4 q+ U Ub 

u *^ y -0.16 

0.45 

60l^ u 

3.11ii 

3.311 1 

3.95 

9 o' 7 +£>->i ) 3 

-1.90 

1133+16 

0 + 2 
°-2 

(3311) 

2711 

q 77+ 0 - 13 

u -' ' -0.15 

U .°0 0.09 

11.11 46 6 6 

2 . 81+2 

9 1 +0.5 
z - i -0.4 

4.25 

n 47+ 0 - 71 
-0.24 


(e.g. iDeshoande fc Rankiri 11999 : lOil fo Sendvki . hoQdtl . 
According to RS75, P 3 = -/VP 3 , where N is the num¬ 
ber of sparks contributing to the drifting subpulse phe¬ 
nomenon observed in a given pulsar and P 3 is the primary 
drift periodicity (distance between the observed subpulse 
drift bands). On the other hand N ~ 2 r Kr p /2h (GSOO). 
Thus, one can write the shielding factor in the form 
77 ~ (1/27t)(P/P 3 ) , which depends only on an easy-to- 
measure primary drift periodicity. Apparently, the shield¬ 
ing parameter 77 should be much smaller than unity. 

The X-ray thermal luminosity is L x = aT^Trr p = 
1.2 x 10 32 (P-i 5 /P 3 )(r]h/r p ) 2 erg/s, which can be com¬ 
pared with the spin-down power E = Iflfl = 3.95/45 x 
10 31 P_i 5 /P 3 erg/s, where I = / 4 5 10 45 g cm 2 is the neu¬ 
tron star moment of inertia (bellow we assume that 
1 .45 = 1). Using equation 0 we can derive the ther¬ 
mal X-ray luminosity and its efficiency as L x = 2.5 x 
10 31 (P_i 5 /P 3 )(P/P 3 ) 2 , or in the simpler form represent¬ 
ing the efficiency with respect to the spin-down power 


h = 0.63 
E 



(3) 


which is useful for comparison with observations. One 
should realize that this equation holds only for thermal 
X-rays from hot spot and cannot be applied neither to 
cooler radiation from the entire stellar surface nor to the 
magnetospheric component. 

We can see that L x in these equations depends only 
on radio observables. It is particularly interesting and im¬ 
portant that this equations does not depend explicitly on 
77 , b and Z, but only on their combination expressed by 
Eq.0. Although one has to be careful whether all our 
assumptions are satisfied in real pulsars, it seems that we 
have found useful, relatively easily testable relationship 
between the properties of drifting subpulses observed in 
radio band and the characteristics of thermal X-ray emis¬ 
sion from the polar cap heated by sparks associated with 
these subpulses. 

Using Eq. 0 we can write the polar cap temperature 
in the form 


T s = (5.1 x 10 6 K)Z 1 / 4 p 1 / 1 5 P " 1/2 



(4) 


where the enhancement coefficient b = B s /B d ~ A pc /Abd, 
A pc = 7 T r pc and Abd = A p is the actual emitting surface 
area (bolometric). Since Aboi can be determined from the 
black-body fit to the spectrum of the observed hot-spot 


thermal X-ray emission, the above equation, similarly to 
Eq.0, depends only on combined radio and X-ray data. 

Although Eq. 0 seems generic in any IAR model in 
which drift rate and the heating rate are determined by the 
same electric field (including pure vacuum gap of RS75), 
Eqs. 0 and 0 considered together describe exclusively 
PSG model. The observational tests of this model are dis¬ 
cussed in section 3 and summarized in Table 1. For PSRs 
B0943+10 and B1133+16, which are the only two pulsars 
in which both P 3 as well as L x and Abd are measured, the 
above equation seems to hold quite well. Interestingly, in 
both cases L x /E ~ 10 -3 and T s ~ 3 MK. 


3. Comparison with observational data 

Table 1 presents the data for two pulsars, which we be¬ 
lieve show clear evidence of thermal X-ray emission from 
the polar caps as well as they have known values of ter¬ 
tiary subpulse drift periodicity. The predicted value of P3 
and/or L x were computed from Eq.0, while the pre¬ 
dicted values of T s were computed from Eq. 0, with 
b = Ape/Aboi determined observationally. 

PSR B0943+10. This is the best studied drifting 
subpulse radio pulsars with P = 1.09 s, P-15 = 3.52, 
E = 10 32 e rg s" 1 , P 3 = 1.86P , P 3 = 37.4P and N = 
P3/P3 = 20 llDeshpande fc Rankinl / l999i'l . It was observed 
by ZSP05, who obtained an acceptable thermal BB fit 
with bolometric luminosity L x = (515'g) x 10 28 erg s _1 
and thus L x /E = (0.491q °1) x 10 -3 . The bolometric po¬ 
lar cap surface area Abd = 10 7 [T s /(3 x 10 6 K)] _ 4 cm 2 ~ 
(llol) x 10 7 cm2 much smaller than the conventional 
polar cap area A pc = 6 x 10 8 cm 2 . This all correspond 
to the best fit temperature T s ~ 3.1 x 10 6 K (see Fig. 1 
in ZSP05). The predicted value of L x /E calculated from 
Eq. 0 agrees very well with the observational data. 
The surface temperature T s calculated from Eq. 0 with 
b = A pc /Abd is also in good agreement with the best fit. 
The shielding factor 77 = 0.09. 

PSR B1133+16. This pulsar with P = 1.19 s, 
P_i5 = 3.7, and E = 9 x 10 31 erg s _1 is almost a 
twin of PSR B0943+10. KPG06 observed this pulsar 
with Chandra and found an acceptable BB fit L^/P = 
(0.77±g;i|) x 10- 3 , A b d = (0.5±g;|) x 10 7 cm 2 and 
T s « 2.8 x 10 6 K. These values are also very close 
to those of PSR B0943+10, as should be expected for 
twins. Using equation 0 we can p r edict P 3 /- P = 27^2 
for B1133+16. Interestingly. iNowakowskl j19961) obtained 
fluctuation spectrum for this pulsar with clearly detected 
long period feature corresponding to about 32P. Most re- 
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cently, WES06 found P3/P = 3±2 and long period feature 
corresponding to (33 ± 3 )P in the fluctuation spectrum of 
PSR B 1133+16. The latte r value seem to coincide with 
that of lNowakowskil lll996lf . as well as with our predicted 
range of P3. We therefore claim that this is the actual ter¬ 
tiary periodicity in PSR B1133+16 and show it in paren¬ 
thesis in Table 1. 

It is worth noting that P3/P3 = 33 ± 3 is quite close to 
37.4 measured in the radio twin PSR B0943+10. Note also 
that the number of sparks predicted from our hypothesis 
is N = P3/P3 = (33 ± 3)/(3 ± 2 ) = nig 5 , so it can also 
be close to 20, as in the case of twin PSR B0943+10. The 
shielding factor 77 = 0.05to'o2- 

4. Conclusions and Discussion 

Within the PSG model of IAR developed by GMG03 we 
derived a simple relationship between the X-ray luminos¬ 
ity L x from the polar cap heated by sparks and the ter¬ 
tiary periodicity P 3 of the subpulse drift observed in radio 
band. In PSRs B0943+10 and B1133+16 for which both 
L x and P3 are known, the predicted relationship between 
observational quantities holds well enough. 

Both the heating and the drifting rate depends on 
the gap potential drop. In this paper we made the point 
that there is continuum of cases between pure vacuum 
gap and the space charge limited flow. The original RS75 
model predicts much too high a subpulse drift rate and 
an X-ray luminosity. Other available acceleration mod¬ 
els predict too low a luminosity and the explanation of 
drifting subpulse phenomenon is generally not clear (see 
ZSP05 for more detailed discussion). Approximately, the 
bolometric X-ray luminosity for t he space charge lim¬ 
ited flow is about (10 -4 4- 10 _5 )P iHarding & Muslimovl 
l 2002 lh and for the pure vacuum gap (RS75) is about 
(10 " 1 4- 1 0~ 2 )E (ZSP05), while for the PSG (GMG03) 
is ~ 10(this paper). The latter model also predicts 
right E x B plasma drift rate. Thus, combined radio and 
X-ray data are consistent only with the PSG model, which 
requires very strong (generally non-dipolar) surface mag¬ 
netic fields. Observations of the hot-spot thermal radiation 
almost always indicate bolometric polar cap radius much 
smaller than the canonical value (at least by an order of 
magnitude). Most probably such a significant reduction of 
the polar cap size is caused by the flux conservation of the 
non-dipolar surface magnetic fields connecting with the 
open dipolar magnetic field lines at distances much larger 
than the neutron star radius. Part of this reduction can 
follow from the fact that discrete sparks do not heat an 
entire polar cap area (but presumably more than a half), 
but uncertainties in determination the size of the polar 
cap are likely to account for this effect. 
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